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Abstract 

The recently developed method to solve the Faddeev equations in coor- 
dinate space is used to study the weakly bound halo nucleus \H(k+n+p). 
The long distances are treated carefully to achieve convergence and high 
accuracy. We use several sets of two-body interactions which reproduce 
the deuteron properties and provide the low-energy A-nucleon scattering 
data close to that of two of Nijmegen potentials. We show that the details 
of the potentials are unimportant unless the accuracy of the hypertriton 
binding energy is required to be better than 50 keV. We find that the 
most significant parameter of the A-nucleon interaction, the singlet s-wave 
scattering length, must be within 10% of 1.85 fm, when the A-separation 
energy from the deuteron is about 130 keV. Other details of the A-nucleon 
interaction are less important for the hypertriton structure. The scattering 
length and effective range are computed for scattering of a A-particle on 
a deuteron. The folding model reducing the three-body problem to a two- 
body problem is investigated in this context and found to be inadequate. 
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1 Introduction 



The last decade has witnessed a rapid development following the discovery of 
nuclear halos |T|. The most carefully studied examples are n Li and 6 He 0. The 
characteristic features of halos are very small particle separation energies and 
very large spatial extension || §]. Simple few-body models are used to describe 
these features and the details of the interactions are not essential even in accurate 
computations of halo properties f| || . 

It is interesting to extend the studies of "ordinary" nuclear halos to systems 
with non-vanishing strangeness. The halo perspective should then be kept and 
the general properties of halos extracted. The lightest nucleus with strangeness 
is the hypertriton (\H). It consists of a neutron, a proton and a A-particle. The 
system may roughly be described as a deuteron and a A-particle bound in a state 
with the binding energy B\ = 0.13 ± 0.05 MeV |7], ||. The deuteron binding 
energy is Bd = 2.224575 MeV and the total three-body binding energy is then 
B = 2.35 MeV. The spin and the parity JJ = 1 + and = \ + , respectively of 
the deuteron and the A-particle, are combined with the corresponding quantum 
numbers of the relative state into the total J 71 " = | + ||. If the deuteron remains 
undisturbed within the hypertriton the total orbital angular momentum can only 
take the values L = or 2. In any case, relative s-states are probably the all- 
dominating components. 

The hypertriton or the approximate A-deuteron two-body system is very 
weakly bound and qualify as the simplest nuclear halo system |l] with finite 
strangeness. The spatial extension of the system is easily estimated by the root 
mean square radius < r 2 > x / 2 ?a h/y/AfiB\ = 10.2 fm, where the reduced mass 
H = mpjrid/ (niA + rrid) = 0.760rriN is about 3/4 of the nucleon mass m^c 2 = 939 
MeV. Thus the A-deuteron radius is 5 times larger than the root mean square 
radius Rd = 1.971 fm of the deuteron. The system appears to be very simple, 
but it is close to the A-dripline and therefore very sensitive to the strength of the 
A-nucleon interaction, which in this way therefore is constrained by the measured 
hypertriton properties. 

The hypertriton offers a mixture of simplicity, sensitivity and fundamental 
interest in connection with the interactions between nucleons and strange parti- 
cles. This combination have attracted considerable (theoretical) attention over 
the last three decades. The interest is reflected in a number of investigations 
where different techniques and models are used, e.g. variational models JR], |TT|j , 
the hyperspherical method 0, [y§, Faddeev calculations JM , 15 1 and simpler 



A-deuteron two-body models [TfJ. Virtually all available interactions have been 
used in these investigations including very sophisticated microscopic interactions. 
Unfortunately the results also exhibit considerable variation for example from the 
(observed) bound three-body structure to the (wrong) unbound A-deuteron sys- 
tem g5|. 

It is also very tempting to use the hypertriton to fine-tune the A-nucleon 



2 



two-body interaction. However, this certainly requires very accurate knowledge 
of both the interactions between all contributing partial waves of the subsystem 
and all contributing couplings to other systems for example like the S-nucleon- 
nucleon system. In addition the rather delicate numerical three-body problem 
also must be mastered to a high accuracy. 

A convenient method to solve the quantum mechanical three-body problem 
was recently formulated ]17| and tested on various systems for example on the 



numerically difficult problem of the occurrence of Efimov states in nuclei [18]. 
It is well-known |5|, § that one of the difficult three-body problems is to give 
an accurate description of a system where one binary subsystem has a loosely 
bound state and the third particle is weakly bound to the other two particles. 
The characteristic size of the A-deuteron system and the related small binding 
energy requires a rather accurate treatment of the wave function to large distances 
up to about 50 fm ||. This is only barely achieved in the best of the previous 
investigations. 

Another interesting observation in this connection is that the large extension 
of the system, which implies high sensitivity to the two-body interactions, at the 
same time provides the simplification that only low-energy properties of the A- 
nucleon interaction is really needed in the description 0, |6| . These properties may 
be expressed in terms of scattering lengths and effective ranges of various relative 
states Thus extraction of other very detailed properties of the interaction 
from the properties of the hypertriton might turn out to be very difficult. 

The purpose of this paper is two-fold. First we want to develop, refine and 



test the new method to solve the three-body problem [17]. An accurate determi- 



nation of the hypertriton structure, including its binding energy, is a challenging 
problem due to the bound subsystem and the weak binding energy 0. The 
present application is an investigation of a "strange" nuclear halo system. In 
this connection it is worth noting that the previous elaborate and complicated 
Faddeev calculations were obtained by using different methods for example by 
working in momentum space. Results from the present new method is desirable 
to get an independent comparison. In addition to the bound state investigation 
we use the occasion to apply the method in a strict three-body computation 
of scattering of a A-particle on a deuteron which subsequently is compared to 
two-body scattering calculations. Secondly, we want to investigate whether the 
low-energy properties of the A-nucleon interaction suffice for the desired accurate 
description of the hypertriton structure. In this connection we need to study the 
effects of the nucleon-nucleon interaction, as reflected in the deuteron properties, 
on the hypertriton structure. 

This paper may be considered as part of a series discussing the general prop- 
erties of three-body halo systems. They all deal with weakly bound and spatially 
extended three-body systems compared to the length and energy of the two-body 
interactions. In the first of these papers we discussed Borromean systems 
where no subsystem is bound H. In the second M we extended the discussion 
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to general three-body systems especially those where two and three-body asymp- 
totics are mixed. In the third ]19| we discussed the effects of the finite spins of the 
particles. In the fourth |2C] we discuss fragmentation interactions of halo nuclei 
in the sudden approximation and especially the effects of final state interactions. 

In this paper dealing with a strange halo, we give in section 2 after the intro- 
duction, a brief sketch of the method. In section 3 we introduce and parametrize 
the all-dominating set of input parameters, i.e. the two-body potentials. In sec- 
tion 4 we describe and discuss the numerical results and in particular the strict 
three-body scattering computations compared to the two-body approximation 
to this three-body problem. Finally, in section 5 we offer a summary and the 
conclusions. 



2 Method 

The method is describe in previous publications, see for example ]5], ^, |19| . There- 
fore we shall here only give a very brief sketch, which should be sufficient to define 
the quantities of interest and the notation used in the following sections. We con- 
sider a system of three particles labeled by % = 1, 2, 3 with corresponding masses 
rrii, coordinates and momenta p^. They interact via the two-body potentials 
Vij. We shall use the Jacobi coordinates, which apart from mass factors, are 
defined as the relative coordinates between two of the particles (x) and between 
their center of mass and the third particle (y). The explicit definitions and the 
corresponding three sets of hyperspherical coordinates (p,a,Q x , Q y ) can be found 
in [||, |^, If. Here p (= y/x 2 + y 2 ) is the generalized radial coordinate and a, in 
the interval [0, 7r/2], defines the relative size of x and y, where Q x and Q y are 
the angles describing the directions of x and y. The volume element is given by 
p 5 dQdp where dfl = sin 2 a cos 2 adadQ x dQ y . 

The total wavefunction \1/ of the three-body system is given as a sum of three 
components ip^) each expressed in terms of one of the three different sets of Jacobi 
coordinates: 

vf^^Wj. (1) 

i=i 

These wavef unctions satisfy the three Faddeev equations 

(T - E)^ + V jk {^ + + ^ k) ) = 0, (2) 

where E is the total energy, T is the kinetic energy operator and {i,j,k} is a 
cyclic permutation of {1, 2, 3}. Each component ip^ is now for each p expanded 
in a complete set of generalized angular functions f^). 

1>® = 4j2 y Ef»(p)* { fM)> ( 3 ) 



4 



where the phase space factor p~ 5 / 2 is extracted. The angular functions are now 
chosen for each p as the eigenfunctions of the angular part of the Faddeev equa- 
tions 

|^A 2 $w + V jk (*$> + *<?> + ^) = ^^A n (p)$«, (4) 

where we choose m = m N and A 2 is the p-independent part of the kinetic energy 
operator defined by 

Explicitly the generalized angular momentum operator A 2 is given by 



-2 "2 

1 <9 2 11 

A 2 = 7— sin a cos a H 1 1 y - 4 (6) 

sin a cos a oa 2 sin a cos 2 a 



~2 -2 

in terms of the orbital angular momentum operators 1^, and l y related to the x 
and y degrees of freedom. 

The radial expansion coefficients f n (p) are solutions to the coupled set of 
differential equations 



2mE 



+ J (Xn(p) + ^))/n + E (- 2P ™'Jp~ - Q»n?J fn> = 0, (7) 



V dp 2 h 2 
where the functions P and Q are defined by 



Pnn'(p) = E / dm®^p,n)^*$(p,n), (8) 



i ,3=1 ' 



Q nn ,( P ) = I dn^*(p,n)^^)(p,n). (9) 
ij=i J °p 

The asymptotic large-distance behavior of the radial potentials can be calculated 
in our case when one of the two-body subsystems has a bound state with the bind- 
ing energy E>d and the corresponding wavefunction 0(r). The result is obtained 
from the behavior of the individual quantities 

Ai _ 4 _W +(0|1 + 3r 9 +r2 |_ Wi (w) 

fti—^ + ^i + s^ + ^W, (ii) 

where the index 1 refers to the lowest A-value. 
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Since the diagonal terms of P vanish, i.e. P nn = 0, the diagonal effective 
radial potential corresponding to the lowest A is then given by 

1 /, 15. ^ 2mB d . o N . . 

-(A: + — ) - Qn - ^ + 0(p- 3 ) . (f 2) 

The diagonal effective potential for higher lying A-values, (A n + 15/4)/p 2 — Q nn , 
instead approaches zero for p — > oo at least as fast as 1/p 3 , see || [l^] . The radial 
couplings arise from the non-diagonal parts of P and Q. They also approach zero 
for p — > oo at least as fast as 1/p 3 . 



3 Two-body potentials 

The present three-body calculation explores how far it is possible to understand 
the hypertriton as a three-body system from its constituent particles and their 
mutual (phenomenological) interactions. The low-energy properties of the in- 
teractions and the interactions at distances much larger than the particle sizes 
are decisive due to the relatively weak binding of the system. We assume that 
the particles maintain their identity and their intrinsic structure only enter via 
their spins and parities. The requirements for the phenomenological two-body 
interactions are then first that they accurately reproduce the low-energy proper- 
ties of the two-body subsystems. Secondly, the interactions must simultaneously 
be simple enough to allow a flexible and accurate treatment of the three-body 
problem. The latter requirement will be met by proper parametrization. 



3.1 The nucleon-nucleon interaction 

The main component of the relevant part of the nucleon-nucleon force corresponds 
to the quantum numbers of the deuteron, i.e. the triplet s and d-states with the 
total angular momentum and parity J 71 = 1 + . However, in the three-body system 
also two-body relative p-states are possible, since corresponding relative p-states 
of the last particle can couple to give the proper total angular momentum and 
parity. Higher orbital angular momenta are negligible. The total isospin of the 
hypertriton is zero and the nucleon-nucleon relative state then also must have 
vanishing isospin. 

We parametrize the isospin zero part of the nucleon-nucleon interaction with 
central, spin-spin, tensor and spin-orbit terms as 

V N n = V c (r) + V ss (r)s m ■ s N2 + V T (r)S 12 + V so (r)l NN ■ s NN , (13) 

where Stvi and Stvi are the spins of the nucleons, s NN = s N i + s N2 , Inn is the 
relative orbital angular momentum and Si 2 is the usual tensor operator. It is 
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sometimes convenient to use a "central" potential including the spin-spin part, 



i.e. 



V {2s + i) = K(r) + Ks ( r )( s | SjV i • s N2 \s) , (14) 

which then depends on the total spin s of the two-body system. To test the 
sensitivity of our results we want to vary the radial shapes of the interactions 
while keeping the deuteron properties and the low-energy scattering properties. 
We choose the radial shapes, V^(r), V^ 3 '(r), Vr(r), V so (r), of each of the terms 
to be either a gaussian Vexp(—r 2 /b 2 ), an exponential Vexp(— r/b) or a Yukawa 
function \^6r _1 exp(— r/b). 

The low-energy scattering data are available as scattering lengths a and effec- 
tive ranges r e of the different relative states. We also have detailed experimental 
information about the deuteron ||21|| , i.e. binding energy E>d = 2.224575(9) MeV, 



root mean square radius Rd =< r 2 > 1 I 2 = 1.971(6) fm, admixture of the d-state 
Pd ~ 4% — 6%, the asymptotic ratio of the d to s-wave component r\d = 0.0256(4) 
and electric quadrupole moment Qd = 0.2859(3) fm 2 . 

In the hypertriton wave function the deuteron component is essential and re- 
producing the deuteron properties is therefore expected to be more important 
than reproducing scattering data. As the s-wave by far is the most dominating 
component we shall first explore this question by using only the central part of 
the interaction and the corresponding deuteron wave function without d-state 
admixture. We then adjust the range and strength of the radial two-body po- 
tential to reproduce the measured values of B d and Rd- We also varied Rd to 
see the influence of the deuteron size on the hypertriton properties. The result- 
ing parameters for these interactions labeled 1 and 2 can be found in the upper 
part of table [l] for gaussian, and exponential radial shapes. The many digits 
in the parameters are necessary to obtain sufficient accuracy on the computed 
quantities. 

We now also include the d-wave admixture in the deuteron. For each radial 
shape we have then four ranges and four strengths to adjust to the four singlet 



and triplet s-state low energy parameters p2[ , binding energy, root mean square 
radius, d-wave admixture, quadrupole momentum and the asymptotic d- to s- 
wave ratio. The parameters for these interactions resulting from this consistent 
treatment are labeled with an additional C and given in the middle part of table [TJ 
In addition we constructed two more sets of potentials which provide slightly 
different deuteron properties in order to check the sensitivity of the resulting 
hypertriton structure to these differences. Those potentials are listed in the lower 
part of table |l|. 

The various properties of the deuteron for each of these nucleon-nucleon in- 
teractions are shown in table ^]. Even for the simplest central interactions we 
obtain scattering lengths and effective ranges for the triplet state in rather close 
agreement with the measured values. The more complicated potentials from 
the middle part of table |I] also reproduce these scattering data rather well. The 
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quadrupole moment, the root mean square radius and the asymptotic d to s-wave 
ratio are reproduced as well although with minor deviations. These quantities are 
nicely determined and some of the remaining small deviations are in fact rather 
difficult to remove due to strong model correlations between the different data. 
The lower part of table ^ shows the properties of the additional sets of potentials. 
We shall use these slightly different properties of the deuteron to study various 
specific correlation with the hypertriton structure. 



3.2 The A-nucleon interaction 

The A-nucleon relative wave function within the hypertriton is dominated by 
s-wave components. In the computations we shall also allow p and d-state con- 
tributions, but we shall neglect the contributions from the higher orbital angular 
momenta. The parametrization of the interaction are again given by the form 

= V} l) (r) + V ss (r)s A ■ s N + V T (r)S 12 + V so (r)l AN • *an , (15) 

where sa, s n , and respectively are the spin of the A-particle, spin of the 
nucleon and the A-nucleon relative orbital angular momentum. The total spin 
is satv = s a + s n and I = 0,1. The central part may have an orbital angular 
momentum dependence indicated by the superscript I. This is to simplify simul- 
taneous reproduction of all the different p-wave scattering lengths. We shall also 
use a parametrization without explicit /-dependence. 

The low-energy properties of the interaction is not directly available from scat- 
tering experiments. However, it is known that the A-nucleon interaction does not 
bind the two-body system and the binding energy of the hypertriton can therefore 
be used to constrain models of this interaction. The role of the hypertriton and 
the A-nucleon interaction for nuclear systems with finite strangeness is in this 
way analogous to the role of the deuteron and the corresponding nucleon-nucleon 
interaction for ordinary nuclei. 

Several models are available for the A-nucleon interaction. We have selected 
two different models developed by the Nijmegen group. They are one-boson 
exchange potentials obtained by using SU(3) symmetry and simultaneous fits to 
all nucleon-nucleon and A-nucleon data EM 1251. We shall here focus on the 



models called Nijmegen F and Nijmegen SC [p3 , [24[] . The computed scattering 
lengths and effective ranges for s-states and the scattering lengths for the p-states 
are shown in table |3[ Using the triplet p-state, this provides 7 constraints on our 
parametrization. For convenience we do not use the singlet p-state in the fits, 
since it is determined (as the s-waves) by the central and spin-spin parts of the 
interaction and the corresponding contributions to the hypertriton wave function 
are small. The low-energy parameters of the Nijmegen models F and SC are 
reproduced by the interactions labeled 1, 3 and 5, 6 respectively. 
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As for the nucleon-nucleon case we determine strengths and ranges of the 
potentials by adjusting to the model values of the s-wave scattering lengths and 
effective ranges and triplet p-wave scattering lengths given at the top of table [| 
The remaining degrees of freedom are constrained by choosing the same ranges for 
some of the terms. These choices are rather arbitrary, but they influence mostly 
the small p-wave part of the A-nucleon relative wavefunction. The resulting 
different sets of interaction parameters (label 1) are collected in table 

These interactions all overbind the hypertriton as we shall see later. The 
binding energy is the decisive quantity determining the structure of the halo and 
we must therefore adjust some of the interaction parameters to reproduce the 
measured B^. This can most efficiently be achieved by changing the scattering 
length of the dominating component in the wave function, i.e. the s-waves. For 
the gaussian potential (Gl) we then simply reduce the strengths of the central 
s-wave potential by 10% resulting in the interactions Glr. For exponential and 
Yukawa shapes we then readjusted the parameters of the potentials to reproduce 
these new slightly altered low-energy scattering data. The resulting parameters 
for these interactions (Ylr and Elr) are also shown in table |j. 

As an alternative choice of the interaction parametrization we do not allow 
now 1-dependence in the central potential leaving 8 parameters and 7 constraints. 
The eighth parameter is constraint by additional fit to the 3 P2 effective range (la- 
bels 3 and 5). We use both the Nijmegen potentials, although mainly Nijmegen 
SC, and reproduce the scattering data with different radial shapes. Again the 
hypertriton binding energy turns out to be too large. We reduce the strength 
of (all partial waves of) the central potential by 10%. The parameters for the 
resulting interactions (labeled G5r, E5r and Y5r) are given in table To study 
the dependence on scattering lengths we varied the singlet and triplet s-wave 
parameters independently. The new sets are labeled "ra" and "rb" and the cor- 
responding interaction parameters for gaussian, exponential and Yukawa shapes 
are also given in table f| (labeled G5ra, G5rb, E5ra, E5rb, Y5ra and Y5rb). 

Finally, we have collected the scattering lengths and effective ranges in table 
|3| for all these A-nucleon interactions. 



4 Numerical results 

The method requires a priori specification of the quantum numbers of the con- 
tributing components in all three Jacobi set of coordinates. A complete basis is 
obtained by including all possible values of these quantum numbers consistent 
with the generalized Pauli principle, parity and the rules for coupling angular 
momenta. The isospin conservation requires that the two-nucleon system is in an 
isospin zero state, i.e. a deuteron-like state. Only very few of the lowest orbital 
angular momenta contribute to the accuracy we need. However, we shall include 
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as many as necessary to obtain the required accuracy. 

In the Jacobi coordinate system, where x refers to the relative coordinate 
between the two nucleons, we need to include orbital angular momenta l x = 
0, 2 with the corresponding spin values s x = 0. The intermediate value l x — 1 
(and the related s x — 1) is not allowed in the two-body system due to parity 
conservation, but the third particle can compensate by being in a relative p- 
state. In both the other two Jacobi coordinate systems, where x refers to the 
relative A-nucleon system, we expect again dominating components with (l x , l y ) = 
(0,0) and admixture of components with (l x ,l y ) = (0, 2), (2, 0), (1, 1). The d- 
wave components can be expected to be similar to the d-wave admixture in the 
deuteron. The resulting set of quantum numbers in our computations are given 
in table [5[ A number of additional channels have also been included to test the 
convergence. Their contributions always turned out to be negligible in the present 
context. 

4.1 Convergence and accuracy 

The large extension or the weak binding of the system relative to a free deuteron 
and A-particle demands an accurate treatment. This means that the wave func- 
tion must be calculated to large distances and small components must be included 
both in the angular Faddeev equations and in the coupled set of radial equations. 
These convergence problems are essential as emphasized previously |12[ [13|] . The 
large distances are handled by analytical extrapolation of the numerically ob- 
tained A-values. This improves significantly both accuracy and speed of the 
computation. 

The separation energy B\ is determined with a relative accuracy of about 10~ 3 
already by including the lowest three values. This is roughly unchanged when the 
interactions result in both somewhat weaker and stronger binding energy. This 
relative accuracy is sufficient for our purpose and in the calculations we therefore 
always included the three lowest A-values in the radial set of equations. 

We first carry out a number of calculations using various combinations of N-N 
and A-N interactions. They are parametrized to fit some of the deuteron prop- 
erties and the low-energy scattering data obtained in the two Nijmegen models. 
The results are shown in table Different radial shapes lead to differences in 
binding energy of up to 340 keV out of the total energy of about 0.31 to 0.65 
MeV. The N-N interactions which best fit the deuteron properties (label C) ex- 
hibit less binding and less deviations resulting from different radial shapes. All 
the interactions considered in the upper part of table |6| produce a significant 
overbinding of the hypertriton. The Nijmegen F model is always closest to the 
measured value. Since the properties depend on the size and binding energy of 
the three-body system, we adjusted the interactions to reproduce the hypertriton 
binding energy. The resulting reduced potentials in table [| are from now on used 
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in the discussion. 

The significance of different components in the wave functions can be seen 
in the middle part of table |6]. The contribution from the p-states amounts to 
« 10 keV compared to the value of B\ « 60-157 keV. One exception is the 
Yukawa interaction in the A-N interaction where the p-states contribute up to 40 
keV. These p-state contributions are not allowed in the isolated deuteron wave 
function, but in the hypertriton they are allowed in combination with the negative 
parity components in the relative A-deuteron wave function. 

The contribution from the d-state in the A-nucleon channel (13 and 14 in 
table ^|) amounts to about 10-25 keV for most interactions and to 55-70 keV for 
the nucleon-nucleon interactions obtained by coupled channel analysis, where the 
tensor interaction effectively is stronger. This d-wave contribution is mediated 
by the d-wave in the deuteron and therefore only present when this admixture 
is included. With d-wave admixture in the deuteron the binding energy of the 
hypertriton is larger first due to the mentioned d-wave contribution in the A-N 
channel and secondly due to the deformation causing larger s-wave contribution. 

The different radial shapes also marginally influence the contributions. The 
wave function moves with decreasing binding energy towards larger distances 
where the effective three-body potential is determined by the scattering lengths 
of the two-body potentials. Thus the detailed radial shapes should then be unim- 
portant. However, different radial shapes of the potentials always add a small 
uncertainty. In the present case of B\ 130 keV and a 2.2 MeV bound binary 
subsystem this uncertainty arises almost entirely from the A-N potential. The 
Yukawa potentials seems to provide up to about 50 keV more binding than the ex- 
ponentials which in turn exceed the binding in gaussian potentials by about other 
50 keV. The Yukawa potentials diverge at short distance which is unphysical and 
rather extreme. The more reasonable finite potentials like gaussians or exponen- 
tials lead to binding energies deviating by less than 40 keV. The deviations due 
to different radial shapes are significantly diminished when the consistent anal- 
ysis of the scattering data are used to constrain the interactions. The resulting 
variations in binding energy then amounts to less than 66 keV between different 
shapes including Yukawa functions. 

The sensitivity against variation of the radial shapes of the potentials is inves- 
tigated by introducing a repulsive core in the central A-N potential. We simply 
use a linear combination of two gaussians where the shortest range corresponds 
to b = 0.53 MeV with a strength of 550 MeV. The longer-range gaussian is then 
used in the usual adjustment procedure where low-energy s-wave scattering data 
are reproduced as described before. We only include central and spin-spin terms. 
The resulting repulsive central potential at r = is then ~ 400 MeV. If the low- 
energy properties of all partial waves should remain unchanged the 1-dependent 
potential for each orbital angular momentum must be individually adjusted. We 
only do this for s and p-waves, see table |j. The effect of this repulsion is of 
course the tendency for the wavefunction to decrease or vanish at short distances 
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between the two particles. 

In an expansion on a basis, as we have used in the present work, this requires 
more basis states with sufficient nodes at short distances. Furthermore, the em- 
phasis is somewhat shifted towards higher angular momenta needed to describe 
the exclusion of the wavefunction from the repulsive region. Both effects increase 
the difficulties in achieving high accuracy and convergence in the numerical com- 
putations. The typical results are a decrease of B\ from 490 keV to 420 keV 
where we already included an energy gain of 22 keV by adjusting the p-wave 
scattering to resemble the result for the case without a repulsive core. For B\ 
around 130 keV the total decrease in binding energy is only 17 keV and the gain 
by p-wave adjustment is reduced to about 17 keV, see table [7]. The uncertainty 
related to a repulsive core is then less than the uncertainty due to different radial 
shapes. 

Finally, we note that the mentioned total uncertainties of around 50 keV are 
extremely small compared to the central values of 10-100 MeV of the two-body 
potentials holding the system together. Thus the dependence on the details of the 
potentials is indeed remarkably small. On the other hand the demand for high 
accuracy at some point necessarily brings in some more details of the potentials. 
In the present case this turns out to be at around 50 keV. Improved accuracy 
can be achieved by an adjustment to the s-wave phase shifts over a range of 
(low) energies. Then the choice of radial shape would be much less important. 
The reason is simply that the decisive low-energy properties deviate less over the 
contributing range of energies. 

4.2 Two-body interactions and the hypertriton structure 

The structure of the hypertriton described as a three-body system is determined 
by the N-N and the A-N interactions. Since the deuteron is a main component 
in this structure the most important parts of the N-N interactions are those 
determining the deuteron structure. The all-decisive property of the deuteron is 
the binding energy and less importantly also the root mean square radius Rd and 
the d-wave admixture. We show in table |5] for different sizes and identical 
binding energies of the deuteron. This A-separation energy decreases about 15 
keV when Rd increases from 1.904 fm to 1.971 fm. The variation of the radial 
shape of the N-N interaction only causes insignificant changes when both Bd and 
Rd are fixed. The inherent uncertainty in these calculations is therefore only a 
few keV, since the N-N interaction and the deuteron properties both are very 
well known. 

The sensitivity to the A-N interaction is higher because the system is close 
to the threshold for binding the A-particle. Then only a little change of the 
potential will change the size and binding dramatically. Since the interaction 
is less known it is of considerable interest to determine the decisive degrees of 
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freedom. This will provide information about which constraints the hypertriton 
properties can place on the interaction. The numerical examples in table |6] show 
that different interactions with the same radial shape and the same s and p low- 
energy scattering data lead to the same hypertriton properties. Different radial 
shapes (i.e., gaussian and exponential) and identical scattering data produce a 
variation of of about 50 keV. 

The same radial shape and a 10% bigger singlet s-wave scattering length (ra) 
increases B\ by about 50-100 keV, respectively for gaussian and exponential 
shapes. A 10% decrease of the triplet s-wave scattering length (rb) decreases the 
Ba only by about 15-30 keV, for both gaussian and exponential shapes of the 
interaction. At least the dependence on the singlet scattering length is significant 
in the present context. Thus, accurate knowledge of the binding energy constrains 
the s-wave scattering lengths and especially the singlet s-wave. The present 
accuracy on B\ allows the claim of about 10% deviation from the value of the 
singlet s-wave scattering length for the interaction Fr. The original two potentials, 
Nijmegen F and Nijmegen SC, certainly both overbind the hypertriton. Their 
singlet s-wave scattering lengths are both too large. 

The separation energy B\ has so far been in focus, but the root mean square 
radius is another important and strongly correlated quantity. It is shown in table 
H for the different interactions. In these tables we also show the mean distance of 
the A-particle from the center of mass of the deuteron under the assumption that 
the deuteron is unperturbed in the hypertriton. The precise definition is given 
by 

2 m N (m A + m d ) , 2 2 

<r Ad>= — — (<P >-2R d ). (16) 

m A m d 

This distance is much larger than the size of the deuteron and the A-particle is 
therefore on average far outside the deuteron. It is easy to verify that the product 
B\- < r\ d > remains approximately constant. The reason for this two-body 
halo dependence is the small separation energy from a "particle" asymptotically 
approaching the deuteron. 

We are now in a position to show the realistic structure of the hypertriton. 
First, we recall that the angular eigenvalues A are closely related to the diagonal 
part of the effective radial potentials, (A n + 15/4)/p 2 — Q nn . The A-spectrum is 
therefore decisive for the spatial structure of the system. An example including 
the five lowest values is shown in fig. 1 for the GC1 and the G5r interactions re- 
spectively for N-N and A-N. The behavior is typical and other realistic potentials 
like the combination GC1 and Y5r or GC1 and E5r would differ by about the 
thickness of the lines already after ~ 3 fm (see fig. la). 

At p = 0, we have the hyperspherical spectrum A = K(K + 4), where K 
is a non-negative even integer. The odd numbers are excluded here for parity 
reasons. The lowest A starts at zero, has a flat minimum and then bends over 
and diverges at large distances parabolically as A = — 2B d m^p 2 jTi 2 . The other 
A-functions reestablish the hyperspherical spectrum asymptotically at large p. 
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Therefore they all approach finite constants at large p. In particular the lowest 
level originating from A = 12 is not far from the lowest A-value at small distance 
close to the attractive pocket. This level, which approaches zero for p —>■ oo, is 
responsible for the amount of relative p-states in the subsystems. 

The radial wave functions corresponding to the three lowest A-values are 
shown in fig. 2. By far the most dominating component is that of the lowest 
A, which at large distances describes a A-particle bound to a deuteron in the 
ground state. The main part of the probability is found between 2 fm and 10 
fm. The other components are very small, but they are essential to achieve the 
accuracy we want. 

A more visual impression of the hypertriton structure is obtained from fig. 3, 
where we plot the probability distribution in the plane defined by the three parti- 
cles. We use the system of coordinates where the center of mass is in the middle of 
the figure, the largest principal moment of inertia is horizontal and the A-particle 
is in the right half. For simplicity, only s-waves are included in the computation 
and we used G2 N-N and G6r A-N interactions (see tables [I] and The prob- 
ability distribution reaches its maximum at the point where the hyperradius is 
equal to 6.20 fm, the hyperangle is equal to 12.27 and the angle between x and 
y is 7r/2. At this point the distance between the nucleons is 1.86 fm and the 
distance between their center of mass and the A-particle is 7.65 fm. 

4.3 Scattering of a A-particle on a deuteron 

The properties of the hypertriton are also accessible through information obtained 
by scattering of A-particles on deuteron. The A-N interaction could be further 
constrained by such scattering measurements. If such information should become 
available and subsequently utilized the corresponding three-body scattering prob- 
lem must be solved. Our method, sketched in section 2, allows computation of 
bound states as well as scattering problems. Both types of computations are 
three-body calculations and in particular the scattering problem does not have 
to rely on two-body or any other approximation. 

The procedure is analogous to the bound state problem. The angular eigen- 
values A are first computed for each value of p. Then the coupled set of radial 
equations are solved with the proper boundary conditions, i.e. exponentially de- 
creasing radial functions f n except for the function corresponding to the lowest 
A-value, which asymptotically describes the deuteron in its ground state. This 
function is at large distances turning into a phase shifted sine-function in the 
variable ky or equivalently kp. When the energy of the A-particle relative to the 
deuteron approaches zero we may in the usual way derive the scattering length 
a p and effective range r p from these phase shifts. 

These quantities are then related to the p degree of freedom as indicated 
by the notation. They must then be transformed to the coordinates describing 
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the relative distance between the deuteron and the A-particle. The resulting 
scattering length and effective range are then given by 



m N (m A + m d ) m N (m A + m d ) 

a A -d = a P \ , r K _ d = r p \ . (17) 

V m A m d V rn A m d 

The hypertriton has angular momentum |, which means that the deuteron 
spin of 1 couples to the spin | of the A-particle. Thus the relative A-deuteron 
state or equivalently the scattering channel has spin |. 

In the present case we aimed directly at computation of the scattering length 
and effective range by using zero energy from the beginning. Then the lowest 
radial wave function is a straight line outside the short-range potential. The 
intersection with the p-axis gives the scattering length a p ( 2 S). The effective 
range r A ^ d ( 2 S) is obtained by direct integration of the expression 



,2 m o m N(m A + m d ) fPm» ,~ 2 2 n 

r A - d { S) = 2\ / (/„ -f ) dp, (li 

V m A m d Jo 



where /o, corresponding to the lowest A- value, is the radial wave function nor- 
malized to unity at the origin. The related wave function f is the radial solution 
without any interaction. It is identical to fo for distances larger than p max or 
equivalently outside the short-range potential. 

Using three A-values as for the bound state calculation, we obtain a scattering 
length accurate up to about 1% and an effective range accurate up to about 
3%. Both quantities are given in the table |6| along with other properties of the 
hypertriton for the various cases considered. The large value of a A ^ d ( 2 S) is a 
natural consequence of the small A-separation energy B A . The connection is 
approximately given by 

h 2 l 

A " 2pa\_ d { 2 S) 1 + r A . d { 2 S)/a A . d { 2 S) ' (19) 

which numerically is fulfilled. The effective ranges are in all cases about 4.0 fm 
± 0.8 fm, i.e. as usual of the same order as the range of the effective potentials. 



4.4 Two-body approximations to the three-body problem 

Reduction of a three-body problem to an effective two-body problem might be a 
computational advantage. It is tempting to try such a reduction for the hyper- 
triton which approximately can be considered to be a A-particle and a deuteron. 
This case appears to be very well suited for a two-body description as also at- 
tempted previously |16| . The first step of the procedure is obviously to construct 



the A-deuteron effective potential by folding the deuteron wave function with the 
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A-N interaction. The remaining steps are then of two-body nature and therefore 
computationally much simpler. 

The interaction between the A-particle and the deuteron is given by 

V Ad (r Ad ) = J dr|0d( r d) (V AN (r A - r N1 ) + V AN (r A - r m )) ip d (r d ) , (20) 

where ip d is the deuteron wave function obtained by solving the Schrodinger 
equation with the nucleon-nucleon potential. The solution ip Ad to the Schrodinger 
equation with this effective two-body potential then gives an approximation to 
the original three-body problem, i.e. the product wave function ^ Ad {Y Ad )%l) d {Y d ). 

The large- distance asymptotic behavior of the correct three-body wave func- 
tion describes a A-particle bound by B A to a deuteron in the ground state. The 
relative coordinate r Ad = y\Jm N {m A + m d ) j '{m A m d ) is then, apart from the mass 
factor, identical to p(~ y) at large distance. After correcting for this scale fac- 
tor on the coordinate, the radial potential in the three-body computation should 
therefore be directly comparable to the two-body potential. 

Fig. 4 shows this comparison between the potential corresponding to the lowest 
A-value and the folding result where for simplicity only s-waves in the A-nucleon 
channel were included. The differences at short distance obviously arise from the 
spatial extension of the deuteron causing the coordinates to differ. The differences 
at large distance are also significant. The folding potential necessarily falls off 
exponentially in accordance with the exponential behavior of the deuteron wave 
function and the short-range A-nucleon potential. On the other hand the three- 
body potential decreases faster than p -3 , but presumably still much slower than 
exponentially. The reason is that the three-body potential also must account 
for polarization effects where the deuteron wave function in the hypertriton is 
distorted away from its ground state shape. Still asymptotically, but not at 
intermediate and short distances, the potentials both correspond to the deuteron 
in its ground state. 

The hypertriton is in this two-body approximation always unbound, but scat- 
tering in the two-body approximation can still be studied and compared with the 
full three-body results. With a binding energy B A pa 130 keV the two-body ap- 
proximation leads to the scattering length and effective range around 4.8 fm and 
3.7 fm respectively. Through eq. fll~9l ) this corresponds to a two-body virtual state 
at around 680 keV, i.e. about 800 keV above the correct energy value. Therefore 
we conclude that the folding model cannot be used to estimate the properties of 
loosely bound systems like the hypertriton. The three-body calculations are then 
necessary to obtain detailed information. 

On the other hand an effective two-body model might sometimes be useful 
in qualitative understanding and in preliminary estimates. We can perhaps get 
some insight by specifying that the folding approximation in eq. fl2"DD basically 
assumes a specific form for the total three-body wave function. The resulting 
energy is therefore only an upper bound and the hypertriton comes out unbound 
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by about 500 keV. Clearly the potential must be more attractive than the folding 
potential if the separation should be reproduced. 

Increasing the depth and/or the range parameter until the measured binding 
energy is reached then provide a two-body model for the hypertriton. The dis- 
tance between A and the deuteron, measured for example as the root mean square 
radius, is then in this model too small due to the too attractive potential needed 
to compensate for the inflexible assumption of the wave function. The smaller 
the A-deuteron separation energy the more extended the system and the better 
is the two-body model. However, processes involving the intrinsic structure with 
one neutron and one proton can not be described. 

An analogous example is the popular n Li described as a three-body system 
consisting of 9 Li surrounded by two neutrons. Again the intrinsic structure can 
not be described in the model. The distances between these three particles are 
of comparable size and detailed and accurate information therefore require three- 
body calculations. Contrary to the hypertriton with one bound binary subsystem, 
the Borromean system 11 Li would be progressively worse described for decreasing 
two-neutron separation energy. 



5 Summary and conclusion 

The general properties of nuclear three-body halo systems have been discussed in 
a recent series of papers. The characteristic features of these bound states are the 
weak binding and the large spatial extension. The hypertriton certainly qualifies 
as a halo system due to the relatively weak interactions between the A-particle 
and the nucleons. Thus strangeness is involved and this system is potentially 
able to provide information of a very different type from that of the usual halo 
nuclei. 

To relate any available experimental information with the basic quantities 
determining the hypertriton structure, it is necessary to perform accurate and 
reliable computations. We solve the coordinate space Faddeev equations by use 
of a recently developed method where the angular part first is obtained for any 
given average distance of the particles. The large distance behavior is known 
analytically and used to improve accuracy and speed of the calculations. The 
coupled radial equations are solved numerically afterwards. The essential large- 
distance behavior is in this way carefully treated as required by the small binding 
and the related large extension. 

We first sketch the method and then we parametrize the crucial input for 
the calculations, i.e. the two-body interactions. We use potentials with central, 
spin-spin, spin-orbit and tensor terms. To study the sensitivity we parametrize 
each of these with gaussian, exponential and Yukawa radial shapes and compare 
the results. For the nucleon-nucleon part we adjust the parameters to constraints 
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obtained from the measured deuteron properties and the s-wave scattering lengths 
and effective ranges. The A-nucleon interaction is much less known. We use 
the low-energy parameters provided by two potentials from the Nijmegen group 
obtained by SU(3) symmetry and fits to all available nucleon-nucleon and A- 
nucleon data. We adjust the parameters to reproduce the corresponding s-wave 
scattering lengths and effective ranges and the p-wave scattering lengths. We 
compare results where d-wave coupling is neglected in the two-body scattering 
analysis used for the parameter extraction. 

The dependence of the hypertriton properties on the various choices are in- 
vestigated in some detail. First the convergence and accuracy of the computa- 
tions are investigated. Three radial equations suffice to an accuracy of 1 part in 
1000. The sensitivity to the nucleon-nucleon interaction is rather low provided the 
deuteron properties are maintained. Both A-nucleon potential models overbind 
the hypertriton by up to half an MeV. Since the binding energy is the crucial 
quantity determining the size of the system, we reduced the attractive strength of 
the A-nucleon central potential by 10%. The hypertriton energy is then roughly 
correct. For a larger binding energy the A-particle is on average closer to the 
deuteron and the hypertriton structure is more sensitive to the details of the 
two-body potentials. Most effects would then be amplified. The sensitivities are 
more realistic when the measured binding energy is used as starting point in the 
investigations. 

The p-wave contribution to the binding energy is small but visible. A typical 
value is 10 keV. It adds a negative parity component to the relative neutron- 
proton wave function. In the three-body system this is compensated by another 
p-wave in the relative A-deuteron wave function. The contribution to the binding 
energy from the d-wave in the A-nucleon relative state is also small but visible 
with a typical value of about 20 keV. 

The most essential part of the A-nucleon interaction corresponds to the singlet 
s-wave. A 10% change of this scattering length changes the binding energy by 
about 70% around the A-separation energy B A 130 keV. Changing towards 
a weaker potential easily then results in an unbound system. For comparison, 
a change of 10% of the triplet s-wave scattering length changes B\ by about 
15 keV. The other scattering lengths are less important in this connection. To 
obtain the correct binding energy B\ we must therefore have a singlet scattering 
length within 10% from the value 1.85 fm. 

The uncertainty arising from the different choices of radial shapes of the A- 
nucleon potential amounts to 50 keV. A repulsive core of up to 550 MeV in- 
troduces an additional inaccuracy of about 15 keV for a A-separation energy of 
130 keV. These uncertainties are significant when compared to the A-separation 
energy of about 130 keV, but extremely small compared to the strengths of the 
two-body potentials. This uncertainty is in fact also very small compared to 2.35 
MeV which is the total hypertriton binding energy including the deuteron part. 
Thus the chosen low-energy scattering data determines the hypertriton structure 
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and its binding energy to an accuracy of about 50 keV. 

The origin of the differences due to the radial shapes evidently must arise 
from finer details of the potentials. The differences would be reduced by adjust- 
ing to phase shifts in an energy range extending beyond that reproduced by the 
scattering lengths and the effective ranges. This corresponds to a better descrip- 
tion of the higher-energy or off-shell behavior. The extreme low-energy scattering 
properties we used in the fits can only predict structures to the accuracy 50 keV 
specified here. The lower the binding energy the better can the properties be 
predicted by the extreme low-energy scattering behavior. 

The three-body scattering problem, where a A-particle is scattered on a 
deuteron, must also contain information about the hypertriton structure or the 
two-body interactions describing this three-body system. Using our method we 
carried out strict three-body calculations of scattering length and effective range. 
The scattering length is closely related to the size of B A and the effective range 
is of the order of the range of the effective three-body radial potential. 

We also tested the two-body approximation obtained by folding the deuteron 
wave function and the A-nucleon interaction. We computed binding energies and 
scattering properties. The approximation is inadequate for the details and the 
accuracy required in the present investigation. 

In conclusion, the hypertriton is a challenging system where a careful treat- 
ment is necessary to connect measured values with the basic interaction involving 
strangeness. The halo structure makes the large distances important and the nu- 
merical work difficult. The details of the interactions enter rather weakly in the 
hypertriton structure which predominantly is determined by a few of the scatter- 
ing lengths. Thus further constraints on the interaction would certainly require 
considerably higher accuracy both in measurements and in calculations. 
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Figure Captions 

Figure 1 The angular eigenvalue spectrum A as function of p corresponding to 
the interactions GC1 and G5r, respectively for the N-N and A-N in- 
teractions. The five lowest A-values are exhibited, a) Short-distance 
behavior of the four lowest A-values for different A-N interaction radial 
shapes. 

Figure 2 The components of the wave function for the lowest three A-values 
corresponding to the calculation in figure 1. 

Figure 3 The probability distribution for the hypertriton in the plane drawn 



through the three particles. The system of coordinates is defined by 
the following constraints: i) the center of mass is in the middle of the 
figure; ii) the principal moment of inertia is along horizontal axes; iii) 
the A-particle is in the right half of the figure. Only s- waves are included 
in the computations. The G2 and G6r interactions, respectively for the 
N-N and the A-N subsystems, are used (see tables |] and |). Regions 
with higher z values are darker. 



Figure 4 The folding potential V\d (dot-dashed curve) compared to the effective 
three-body radial potential V e ff (solid curve) for the interactions GC1 
and G6r. The coordinate used is p and the mass corrected distance be- 



tween the A-particle and the deuteron y = r^J (wAWd)/mjv(wA + md)- 
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Table 1: Parameters of the nucleon- nucleoli (N-N) interaction. Strengths V in 
MeV and ranges r in fm of the singlet and triplet central (or central (c) and spin- 
spin (ss)), tensor (t) and spin-orbit (so) parts are listed. The radial shapes are 
gaussian (G), exponential (E) and Yukawa (Y), respectively. The upper part of 
the table refers to a deuteron in the 3 Si state; the lower part refers to a deuteron 
in the coupled 3 Si- 3 Di two-nucleon state. 
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Table 2: Properties of the deuteron for the various interactions in table [T]. The 
columns contain the root mean square radius < r 2 > 1//2 , the ci-state probability 
Pd, the electric quadrupole momentum Qd, the asymptotic d to s-wave ratio Tja, 
and the s-wave scattering lengths a and effective ranges r e . The deuteron binding 
energy is Bd = 2.224575 MeV. The first row shows the experimental values. Note 
that in our sign convention a negative scattering length corresponds to a bound 
state. 
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Table 3: Scattering lengths and effective ranges in fm for the two selected A- 



nucleon (A-N) models |23| , g4|] . The rows labeled Fr and SCr represent the values 
obtained by a 10% reduction of the s-wave central strength of the gaussian inter- 
actions, which otherwise reproduces the values for these quantities obtained from 
the two models. The rows labeled ra and rb represent variations where either the 
singlet S or the triplet S scattering lengths are changed from the values of Fr and 
SCr by 10%, respectively up or down. Finally, the two last rows represent the 
the s-wave low-energy data of the models SC and SCr respectively, reproduced 
by using gaussian interactions with only central and spin-spin terms, with (label 
cp) and without repulsive core, see table |j. 
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-0.128 


0.061 


-0.551 


Y5 


2.78 


2.88 


1.41 


3.42 


0.432 


-0.059 


0.062 


-0.525 


G5r 


2.27 


3.10 


1.15 


3.36 


0.405 


-0.068 


0.086 


-0.415 


E5r 


2.27 


3.10 


1.15 


3.54 


0.407 


-0.075 


0.085 


-0.451 


Y5r 


2.27 


3.10 


1.16 


3.85 


0.397 


-0.066 


0.085 


-0.494 


G5ra 


2.50 


3.11 


1.15 


3.36 


0.451 


-0.067 


0.085 


-0.416 


E5ra 


2.49 


3.10 


1.15 


3.48 


0.446 


-0.072 


0.086 


-0.442 


Y5ra 


2.49 


3.10 


1.15 


3.57 


0.428 


-0.065 


0.086 


-0.457 


G5rb 


2.27 


3.11 


1.04 


3.36 


0.408 


-0.064 


0.087 


-0.366 


E5rb 


2.27 


3.10 


1.04 


3.64 


0.407 


-0.067 


0.086 


-0.413 


Y5rb 


2.27 


3.13 


1.05 


3.87 


0.399 


-0.053 


0.087 


-0.450 


G6, G6cp 


2.78 


2.88 


1.41 


3.11 


0.434 


-0.230 


-0.230 


-0.230 


G6r, G6rcp 


2.27 


3.10 


1.15 


3.36 


0.407 


-0.201 


-0.201 


-0.201 
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Table 4: Parameters of the A-nucleon interaction (A-N). Strengths V in MeV 
and ranges r in fm of the central (c), spin-spin (ss), tensor (t) and spin-orbit 
(so) parts are listed. The central part V} l '(r) depends on the angular momentum 
/ = 0,1. The radial shapes are gaussian (G), exponential (E) and Yukawa (Y), 
respectively. The corresponding interactions reproduce the scattering lengths of 
the Nijmegen F model |23| (labels 1 and 3), the Nijmegen SC model |24| (labels 
5 and 6), the reduced Nijmegen models (label r) and the modified potentials 
(label ra and rb) corresponding to Fra, Frb, SCa and SCb in table |3], see also 
the text. The interactions with labels 3 and 5 all include d-wave coupling in 
the analysis. Label c refer to interactions with a fictitious repulsive term in the 
central (1=0) part with parameters V cr = 550 MeV and r cr = 0.53 fm. Label 
p refers to interactions reproducing the same sets of low-energy data (including 
p-wave scattering lengths) as the interactions labeled G6 and G6r, respectively 
(see table |3T). 



A-N 


v c 


r (l=0) 
c 


v^ =1 ) 

y c 


(i=i) 

c 


V ss 


'"ss 


v t 


n 


V so 


' so 


Gl 


-24.354 


1.4986 


-5.413 


1.4986 


1.931 


1.5998 


1.884 


1.4986 


-2.347 


1.4986 


Yl 


-50.833 


0.8726 


-10.010 


0.8726 


2.529 


1.0764 


3.296 


0.8726 


-4.109 


0.8726 


G3 


-24.272 


1.4972 


-24.272 


1.4972 


2.210 


1.5486 


2.053 


1.7278 


-3.644 


2.0488 


Glr 


-21.919 


1.4986 


-5.413 


1.4986 


1.931 


1.5998 


1.884 


1.4986 


-2.347 


1.4986 


Elr 


-86.859 


0.5678 


-28.426 


0.5710 


6.199 


0.6503 


7.696 


0.5710 


-14.828 


0.5710 


Ylr 


-44.093 


0.8912 


-10.010 


0.8726 


1.997 


1.2141 


3.296 


0.8726 


-4.109 


0.8726 


G5 


-28.280 


1.361 


-28.280 


1.361 


4.838 


1.8221 


4.364 


1.4401 


-8.445 


1.5902 


E5 


-97.431 


0.5446 


-97.431 


0.5446 


25.420 


0.6049 


15.604 


0.5740 


-4.964 


1.0270 


Y5 


-56.291 


0.8119 


-56.291 


0.8119 


4.684 


1.4193 


2.347 


1.1822 


-2.214 


1.7113 


G5r 


-25.452 


1.361 


-25.452 


1.361 


4.838 


1.8221 


4.364 


1.4401 


-8.445 


1.5902 


E5r 


-98.351 


0.5218 


-98.351 


0.5218 


18.896 


0.6812 


16.122 


0.5532 


-2.565 


1.2230 


Y5r 


-50.654 


0.8145 


-50.654 


0.8145 


9.725 


1.0234 


0.922 


1.5597 


-2.733 


1.6423 


G5ra 


-25.185 


1.3759 


-25.185 


1.3759 


5.122 


1.8748 


3.745 


1.4837 


-9.056 


1.5667 


E5ra 


-98.693 


0.5236 


-98.693 


0.5236 


19.195 


0.7137 


16.780 


0.5465 


-2.302 


1.2634 


Y5ra 


-53.974 


0.7977 


-53.974 


0.7977 


8.172 


1.1573 


0.318 


2.2190 


-5.410 


1.3073 


G5rb 


-25.506 


1.3388 


-25.506 


1.3388 


5.443 


1.8284 


4.618 


1.4062 


-8.935 


1.5383 


E5rb 


-96.507 


0.5172 


-96.507 


0.5172 


22.100 


0.6752 


17.087 


0.5391 


-2.390 


1.2215 


Y5rb 


-51.068 


0.8006 


-51.068 


0.8006 


10.151 


1.0598 


0.444 


1.9396 


-5.425 


1.3061 


G6 


-28.644 


1.3588 


-28.644 


1.3588 


4.694 


1.8178 












G6c 


-152.876 


0.9853 


-152.876 


0.9853 


5.284 


1.7492 












G6cp 


-152.876 


0.9853 


-236.971 


0.8468 


5.284 


1.7492 












G6r 


-25.723 


1.3614 


-25.723 


1.3614 


4.574 


1.8368 












G6rc 


-151.130 


0.9729 


-151.130 


0.9729 


5.094 


1.7750 












G6rcp 


-151.130 


0.9729 


-319.318 


0.7722 


5.094 


1.7750 
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Table 5: The contributing three-body channel quantum numbers for the hyper- 
triton ground state. The upper part of the table refers to the Jacobi coordinates 
with the A-particle as spectator, while the lower part refers to the other two 
Jacobi systems where one of the nucleons is the spectator. The quantity t x is the 
total isospin of the two-particle system connected by the relative coordinate x. 



Channel 


Config. 




ly 


L 


s x 


s 


tx 


1 


3 Si 











1 


1/2 





2 


'Pi 


1 


1 








1/2 





3 


'Pi 


1 


1 


1 





1/2 





4 


3 £>1,2 


2 





2 


1 


3/2 





5 


'So 














1/2 


1/2 


6 


3 ^ 











1 


1/2 


1/2 


7 


'Pi 


1 


1 








1/2 


1/2 


8 


'Pi 


1 


1 


1 





1/2 


1/2 


9 


* Pj 


1 


1 





1 


1/2 


1/2 


10 


3 p. 


1 


1 


1 


1 


1/2 


1/2 


11 


3 p. 


1 


1 


1 


1 


3/2 


1/2 


12 


3P, 


1 


1 


2 


1 


3/2 


1/2 


13 


3 ^ 





2 


2 


1 


3/2 


1/2 


14 


3 £l,2 


2 





2 


1 


3/2 


1/2 
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Table 6: The \H binding energy in MeV relative to the deuteron for different 
two-body interactions. The A-N interactions of the upper part (double separation 
line) reproduce the low-energy data of the F and SC Nijmegen models, while those 
of the lower part refer to the reduced ones (Fr and SCr), see tables [L|-|4]. Labels 
a and /3 indicate which channels of table [5] are included: "all except 13 and 14", 
and "only 1,4,5,6", respectively The f3 combination excludes all p- and d-states 
in the A-N subsystem. No greek labels means that all channels are included. The 
\H and the A-deuteron rms radii, and the A-d 2 S- scattering length (negative in 
our notation) and effective range are also showed. All lengths are in fm. 



N-N 


A-N 


B A 


<r 2 > l/2 


^ T 2 ^1/2 


a A - d ( 2 S) 


r A - d ( 2 S) 


G3 


Gl 


0.31 


4.03 


7.72 


-11.5 


3.2 


G3 


Yl 


0.53 


3.06 


5.52 


-9.0 


2.8 


Y3 


Gl 


0.39 


3.96 


7.48 


-10.8 


3.4 


Y3 


Yl 


0.65 


3.21 


5.75 


-8.9 


3.2 


GC1 


G3 


0.22 


4.53 


8.79 


-13.0 


3.3 


EC1 


G3 


0.23 


4.36 


8.41 


-12.5 


3.2 


GC1 


G5 


0.37 


3.71 


6.95 


-10.5 


3.1 


GC1 


E5 


0.44 


3.44 


6.35 


-9.7 


2.9 


GC1 


Y5 


0.51 


3.26 


5.92 


-9.2 


2.8 


EC1 


G5 


0.40 


3.58 


6.67 


-10.2 


3.0 


Gl 


Glr 


0.064 


7.29 


14.76 


-22.6 


3.1 


Gl 


Glr («) 


0.063 


7.32 


14.81 


-22.7 


3.1 


Gl 


Glr W 


0.055 


7.62 


15.44 


-24.1 


3.1 


Gl 


Ylr 


0.158 


5.15 


10.19 


-14.9 


3.0 


Gl 


Ylr W 


0.157 


5.17 


10.23 


-14.9 


3.0 


Gl 


Ylr W 


0.117 


5.40 


10.72 


-15.2 


3.0 


El 


Elr 


0.114 


5.86 


11.72 


-17.2 


3.0 


El 


Elr ( Q ) 


0.113 


5.89 


11.77 


-17.3 


3.0 


El 


Elr W 


0.101 


6.14 


12.32 


-18.2 


3.0 


G2 


Glr 


0.046 


8.03 


16.28 


-26.1 


3.3 


E2 


Glr 


0.049 


7.87 


15.95 


-25.4 


3.2 


G3 


Glr 


0.091 


6.55 


13.14 


-19.3 


4.8 


G3 


Elr 


0.145 


5.37 


10.66 


-15.7 


4.6 


G3 


Ylr 


0.172 


5.01 


9.89 


-15.6 


4.9 


GC1 


G5r 


0.108 


5.88 


11.73 


-17.8 


3.3 


GC1 


E5r 


0.161 


5.09 


10.01 


-14.9 


3.2 


GC1 


Y5r 


0.174 


4.95 


9.71 


-14.4 


3.2 


GC1 


G5ra 


0.166 


5.06 


9.96 


-14.8 


3.3 


GC1 


E5ra 


0.235 


4.41 


8.53 


-12.7 


3.1 


GC1 


Y5ra 


0.247 


4.33 


8.37 


-12.5 


3.2 


GC1 


G5rb 


0.084 


6.25 


12.52 


-19.3 


3.3 


GC1 


E5rb 


0.130 


5.48 


10.87 


-16.3 


3.2 


GC1 


Y5rb 


0.149 


5.23 


10.32 


-15.3 


3.2 


EC1 


G5ra 


0.176 


4.87 


9.53 


-14.1 


3.3 


EC1 


G5rb 


0.100 


5.96 27 


11.89 


-18.1 


3.2 



Table 7: The hypertriton binding energy in MeV relative to the deuteron for 
different two-body interactions. The N-N interaction is always the one labeled 
GC1, see table [l]. The A-N interactions all reproduce the s-wave low-energy data 
of the Nijmegen models SC and SCr, see table [|. The interactions labeled G6cp 
and G6rcp reproduce the same s-wave low-energy data and p-wave scattering 
lengths as the interactions labeled G6 and G6r. The values with the (a) include 
the energy gain obtained adjusting the p-wave scattering lengths to resemble the 
result for the case without a repulsive core. The first column indicates the partial 
waves included in the computations. 



partial waves 


A-N 


B A 


A-N 


B A 


A-N 


B A 


s 


G6 


0.442 


G6c 


0.303 


G6cp 


0.303 


s,p 


G6 


0.484 


G6c 


0.365 


G6cp 


0.387 


s,p,d 


G6 


0.490 


G6c 


0.396 


G6cp 


0.418 a 


s,p,d,f 


G6 


0.491 


G6c 


0.397 


G6cp 


0.419 a 


s,P,d,f,g 


G6 


0.491 


G6c 


0.397 


G6cp 


0.419 a 


s 


G6r 


0.137 


G6rc 


0.083 


G6rcp 


0.083 


s,p 


G6r 


0.166 


G6rc 


0.119 


G6rcp 


0.136 


s,p,d 


G6r 


0.170 


G6rc 


0.125 


G6rcp 


0.152 a 


s,p,d,f 


G6r 


0.170 


G6rc 


0.126 


G6rcp 


0.153 a 


s,P,d,f,g 


G6r 


0.170 


G6rc 


0.126 


G6rcp 


0.153 a 
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